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ABSTRACT 



We examined the gravitational contraction of isothermal molecular cloud cores with slow 
rotation by means of two-dimensional numerical simulations. Applying a sink-cell method, we 
followed the evolution of the cloud cores up to the stages at which most of the matter accretes 
onto the central region (i.e., a protostar and a rotationally-supported circumstellar disk). We 
show that both an infalling disklikc envelope and a rotationally-supported disk around the 
central star arc natural outcome of the gravitational contraction of a prolate cloud core with 
slow rotation. The early evolution of the infalling envelopes resembles sheet models recently 
proposed by Hartmann and coworkers. In the infalling disklike envelope, the radial profiles of the 
density, radial velocity, and azimuthal velocity can be approximated by p cx r^-^-^, Vr <x r~°'^, 
and oc r~^, respectively. The fate of the infalling envelopes is also discussed. 



Subject headings: accretion, accretion disks — hydrodynamics — ISM: clouds — rotation — 
stars: formation 
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1. Introduction 

Recent high-resolution interferometric observations have suggested that diskhke envelopes are common 
structures around embedded young stars (e.g., Hayashi, Ohashi, & Miyama 1993; Ohashi et al. 1996). In 
the disklike envelopes, infall motions seem to dominate over the rotation around the young stars, indicating 
that the disklike envelopes are not rotationally-supported. Thus, embedded young stars are most likely 
to be undergoing dynamical accretion from their disklike envelopes onto quasi-static protostellar cores 
and rotationally-supported circumstellar disks. The typical infalling disklike envelope has a size of the 
order of 10^ AU, which is an order of magnitude larger in size than a rotationally-supported circumstellar 
disk observed in the dust emission. The infalling envelopes are expected to play a crucial role in the 
formation processes of stars as well as their surrounding rotationally-supported disks. It is thus important 
to understand the formation and evolution of the infalling disklike envelopes. 

Several authors have proposed the formation mechanisms of such infalling disks around young stars 
by means of numerical hydrodynamics. Galli & Shu (1993) examined the gravitational contraction of 
magnetized spherical cloud cores [see also Boss (1987) and Yorke & Bodenheimer (1999) for the evolution 
of non- magnetized rotating spherical cloud cores]. They showed that the magnetic field impedes the 
contraction perpendicular to the field lines and thus the disklike structures develop in the infalling envelopes. 
In their model, the magnetic field is essential in the formation of the infalling disks. Hartmann et al. (1994) 
proposed another process for the formation of the infalling disks. They showed that even in the absence 
of the magnetic field, initially sheetlike cloud cores collapse to form highly disklike envelopes which have 
similar structures to Galli & Shu (1993) 's disks. Their calculations indicate that the cloud geometry plays 
an important role in the formation of the infalling envelopes. 

Radio observations have shown that most molecular cloud cores are elongated, implying on statistical 
grounds that most of them are prolate (Myers et al. 1991). Bonnell, Bate, & Price (1996) thus examined 
the gravitational contraction of non-rotating prolate cloud cores. They found that in a prolate cloud core 
close to virial equilibrium, the pressure force retards the contraction preferentially along its minor axis. As 
a result, an extended disklike structure, qualitatively similar to Hartmann et al. (1994) 's sheet model, is 
formed. 

Although these calculations have demonstrated the formation processes of the infalling disks in 
molecular cloud cores, it is still not clear how the infalling disks form and evolve subsequently and how 
rotationally-supported circumstellar disks form in the infalling disklike envelopes. 
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In this paper, we examine the gravitational coUapse of prolate cloud cores with slow rotation. We show 
that both the infalling disklike envelope and rotationally-supported circumstellar disk are natural outcome 
of the gravitational contraction of a prolate cloud core. In §^ we describe our model and numerical method. 
Numerical results are given in §|[ In §|[ we discuss the evolution of the mass infall rate in the disklike 
envelope. Finally we discuss implications to observations. 



2. Models and Numerical Methods 

Radio observations have revealed that most molecular clouds have more or less filamentary structures. 
Recent numerical simulations have also shown that such filaments are ubiquitous features in self-gravitating 
molecular clouds (e.g., Klessen, Burkert, & Bate 1998). According to a linear stability analysis, the linear 
perturbation grows fastest in time when its longitudinal wavelength is several times as long as the effective 
cloud diameter. It is worth noting that in the absence of the strong magnetic fields and fast rotation, the 
filamentary cloud is unstable to linear perturbations whose wavelengths are more than twice the effective 
cloud diameter (Matsumoto, Nakamura, & Hanawa 1994). Hence, the fragments tend to become prolate at 
the epoch of the fragmentation. We investigate the evolution of such prolate fragments formed through the 
fragmentation of filamentary clouds. 

The model cloud is the same as that of Matsumoto, Hanawa, & Nakamura (1997). We consider an 
isothermal cylindrical cloud in which the density is uniform along the cylinder axis. We assume that the 
cloud is in hydrostatic equilibrium in the radial direction and that the cloud is rotating around the cylinder 
axis. The ratio of the centrifugal force to the pressure force, (3, is assumed to be spatially uniform. The 
radial distributions of the density and velocity are then given by 



„2 



V 



and 



V = {Vr,V^,Vz) = 



r2 \"'/' 



(2) 



respectively, where = (1 + (3)c'^/{47TGpo), — 2TTGpof3/{l + (3), Cg is the isothermal sound speed, G is 
the gravitational constant, and po is the initial central density. It should be noted that in this model, most 
of the matter rotates nearly uniformly {v^ cx r) and only the outer low-density region with a small mass 
rotates differentially {v^p = const.). 

The initial model of our simulation is obtained by adding a velocity perturbation onto the equilibrium 
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model. The velocity perturbation changes sinusoidally with the wavelength of A^, 

5v^ = vosm{2nz/X^) , (3) 

where vq is the amplitude of the perturbation and is set to 5% of the sound speed. In the following, we 
take po and Cg as the units of density and velocity, respectively. Then, the initial model is specified by the 
two parameters of (3 and Xz- For the typical models, the wavelength Xz is taken to be 15.5 Cs/ \/2TrGpQ, 
which is almost equal to those of the fastest-growing linear modes for non-rotating or slowly-rotating clouds 
(Matsumoto et al. 1994). Observationally, most molecular cloud cores are likely to be rotating quite slowly 
(Goodman et al. 1993). Thus, we adopt the small values of /3 (= — 0.1). In Table ^ we summarize the 
initial conditions of the models calculated in this paper. 

We solved numerically the time evolution of the model cloud with the spatially and temporally 
second-order- accurate hydrodynamical code described by Matsumoto et al. (1997). Since our model keeps 
mirror symmetry at z = and Xz/2 throughout the evolution, we restricted the computational domain 
in the interval of < z < Xz/2. We set the fixed boundary conditions at r = 0, i?max and z = 0, Zj^ax, 
where i?max and .^max denote the lengths of the computational box in the r- and z-directions, respectively. 
For all the models, we take i?niax = S^max = ^z- The effect of the fixed boundary is very small. We have 
calculated the evolution of the models with larger i?max and have confirmed that the numerical results are 
not changed. 

Many numerical simulations have shown that the gravitational collapse of an isothermal cloud 
accelerates in the central high-density region. Thus, the characteristic length and time scales shorten as the 
collapse proceeds, this phase being referred to as runaway collapse phase. Therefore, it is difficult to resolve 
the central high-density region sufficiently at the late stages of the computation even if the cloud evolution 
is followed with very fine grids. In real evolution, as the collapse progresses, a protostellar core is formed at 
the center and thereafter the envelope gas begins to accrete onto the central core, this phase being referred 
to as the accretion phase. 

To pursue the evolution of the infalling envelope gas, we adopted a sink-cell method when the central 
density reaches a reference density psink [see, e.g.. Boss & Black (1982) and Tomisaka (1996)]. In this 
paper, we take Psink = 10** po and the innermost region oi < r < 5Ar and < ^ < 5Az is assigned to the 
sink cells, where Ar and Az denote the grid spacings in the r- and z-directions, respectively. The excess 
density (p — Psink) and the corresponding momentum are removed from the sink cells in every time step. 
The removed mass is added to a central point mass that affects the matter in the envelope only through 



- 6 - 



gravitation. 

As mentioned above, we take po = 1, Cs = I, and 2TrG ~ 1. Then, the units of time, length, and mass 
are given as 

(„ = (2»Gpo)-» = 

ru = c,(27rGpo)""' 



= l.lxlO^AU — r L^. J , (5) 

V 0.2km s-i/ \106cm-3/ ^ ^ 



"lu = Perl 



^ V0.2kms-i/ Vl06cm-3/ ^ ^ 

respectively, where tiq = Po/ifJ-m-n) is the number density, /i is the mean molecular weight and taken to be 
2.3, mfj is the mass of a hydrogen atom. 



3. Numerical Results 

3.1. Prolate Cloud Cores without Rotation 

In this subsection, we explore the evolution of non-rotating models. In the following, we mainly 
show the numerical results of model B as a typical example. In this model, the wavelength Az is almost 
equal to that of the fastest growing linear perturbation (e.g., Nakamura, Hanawa, & Nakano 1993). This 
model is essentially the same as model D of Nakamura, Hanawa, & Nakano (1995) who followed the cloud 
contraction during the runaway collapse phase. We extend their calculations to considerably later stages, 
i.e., the accretion phase, by applying the sink-cell method. We pursued the accretion phase evolution until 
most of the envelope gas accretes onto the central point mass. 

Figure |l] shows the density and velocity distributions in the r — z plane at eight different stages. The 
first three panels show the cross sections of the cloud during the runaway collapse phase, while the others 
show those during the accretion phase. At the early stages, the cloud collapses preferentially towards the 
z = plane (Figs. |l]b and|l|c). As the collapse proceeds, the central high-density region becomes spherical, 
whereas the low-density region remains prolate. The central spherical region continues to contract towards 
the center. Such evolution is the same as that of model D by Nakamura et al. (1995). 

When the central density reaches W^po, the sink cell method is applied and the central point-like 
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object forms. The mass of the central point-hke object monotonously increases with time. Since the vertical 
flow is dominant near the center, the disklike envelope is formed around the central point-like object and 
grows in radial extent (Figs. [|d and |l|e). A shock wave is also formed at the disk surface. By the stage at 
which the mass of the central point- like object reaches about 0.7Mt, the infalling disk extends radially up 
to r ~ 4 — 5ru, where Mt denotes the total cloud mass. In the disk, the gravitational force dominates the 
pressure force in the radial direction, while in the vertical direction, the pressure force becomes comparable 
to the gravitational force. Therefore, the gas flow progresses mainly in the radial direction inside the disk. 
Such an non-equilibrium disk resembles pseudo-disks found by Galli & Shu (1993), Hartmann et al. (1994), 
and Bonnell ct al. (1996). 

After the mass of the central point-like object reaches about 0.7Mt, the radius of the infalling disk 
becomes shorter because of dominant radial infall (Figs. |l]f and |^g). The shock front lifts up vertically 
because of dominant pressure force in the vertical direction. Thereafter, the shock front reaches the z-axis 
and thus the matter in the infalling disk collapses to form a thin needle-like (prolate) object along the 
z-axis (Figs. |l|h and|l|i). 

We also followed the evolution of models A and C to examine the effects of A on the formation and 
evolution of the disklike envelope. It is found that the evolution is qualitatively similar to that of model B. 
In both models, the disklike envelopes are formed and then they evolve into thin needle-like objects. For 
the model with larger A, the vertical infall is more dominant and thus the disk is more extended in the 
radial direction. 



3.2. Prolate Cloud Cores with Slov^r Rotation 

In this subsection, we explore the evolution of slowly-rotating clouds. In the following, we mainly show 
the numerical results of model E as a typical example. In this model, the wavelength A^ is almost equal to 
that of the fastest growing linear perturbation. As shown below, the evolution of model E is qualitatively 
similar to that of model B except near the center, i.e., the cloud collapses to form an infalling disk in which 
the gas flows mainly in the radial direction. 

Figure ^ is the same as Figure |^ but for model E. During the runaway collapse phase, the central 
high-density region becomes slightly oblate because of rotation, whereas the low-density region remains 
prolate. In the central region, the rotation is not important in the cloud support. Therefore, the central 
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region continues to contract. 

After the central core is formed, the diskhke envelope is formed and then expands in radius (Figs. |^ 
and|^). A shock wave is also formed at the disk surface and propagates towards the z-axis. Such evolution 
is essentially the same as that of model B. In the disk, the rotation is not important and therefore the gas 
flow progresses primarily in the radial direction. As the contraction proceeds, the centrifugal force becomes 
important at the innermost region and a rotationally-supported disk is formed there. 

After the mass of the central region reaches about 0.7Mt, the radius of the infalling disk becomes 
shorter owing to dominant radial infall (Figs. ||f and ^). Then, the matter in the infalling disk collapses 
to form a funnel in which the centrifugal force balances with the gravitational force in the radial direction 
(Figs. ^1 and ^). In Figure ^, we show the enlargements of the central region at the same stages as the 
last three panels of Figure |[ A new shock front is also formed there and propagates in the radial direction 
as the gas with higher specific angular momentum accretes. Such a structure resembles a bipolar cavity 
which is often thought to be formed by the outflow from the central young star. This structure might play 
a significant role in the guidance of outflow from the central star at the late stages of the evolution. 

The rotationally-supported disk is encased in two accretion shock fronts, both of which are located 
several scale heights above the equatorial plane (Figs, and^). Such structure seems to be similar to that 
of Yorkc & Bodenheimer (1999). Shock formation is explained as follows. In the infalling disk, the accretion 
flow is dominant in the radial direction. The radial infall stops at the funnel owing to the centrifugal barrier 
and thus the outer shock front is formed. Thereafter, the matter in the rotationally-supported disk collapses 
in the vertical direction to achieve hydrostatic equilibrium. This infall produces the inner shock front. In 
the innermost disk, angular momentum transfer is expected to be important in the accretion flow onto the 
central core, although it is not included in the present calculations (see Yorke & Bodenheimer 1999). 

We also pursued the evolution of models D and F to examine the effects of f3. It is found that the 
evolution is qualitatively similar to that of model E, although the size of the rotationally-supported disk 
depends on /?, i.e., when (3 is larger, the rotatinally-supported disk is wider. 

3.3. Effects of Cloud Geometry (Prolate vs. Oblate Cloud Cores) 



As shown in §3.1 and §3.2, a prolate cloud core collapses to form an infalling disk in which the accreting 
gas flows mainly in the radial direction. In this subsection, to clarify the physics of the formation process 
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of the infalling disk, we calculate the evolution of two representative models: (1) a prolate cloud and (2) an 
oblate cloud. In the following, we consider non-rotating prolate and oblate clouds which have Gaussian-like 
density profiles as 



Po exp 



Rl Zl 



(7) 



At the initial state, the cloud is assumed to be static for simplicity. This model has the two parameters of 
i?o and Zq. 

Figures ^ and ^ show the evolutions of the density distributions in the r — z plane for the 
prolate and oblate clouds, respectively. The prolate and oblate clouds have the model parameters of 
{Rq,Zq) = (1.66ru, 2.49ru) and (2.40ru, 1.60ru), respectively. In both models, the ratio of the major to 
minor axes is set to 1.5 and the ratio of the thermal energy to the gravitational energy is equal to a ^ 0.33. 
The total cloud masses are equal to Mt — 1.44 x lO^mu and 1.40 x lO^mu for the prolate and oblate 
clouds, respectively. The grid spacing is set to Ar ~ Az — 0.02r„. The number of grid points is taken 
to be Nr 'X Nz = 512 x 1024 In Figures ^ and ||, the first two panels show the cross sections during the 
runaway collapse phase, while the rest are during the accretion phase. In both models, the contraction 
progresses preferentially along the major axis. Therefore, during the accretion phase, the initially prolate 
cloud collapses into an oblate core, whereas the initially oblate cloud collapses into a prolate core. 

This evolution is explained as follows. At the initial state, the pressure force is less effective along the 
major axis because of shallower density gradient. In other words, the net attractive force is stronger along 
the major axis. Therefore, the cloud collapses preferentially along the major axis. 

We also followed the evolution of the models with different axis ratios and a. It is found that as long 
as the cloud is close to virial equilibrium, the evolution is qualitatively similar to those shown here (see also 
Bonnell et al. 1996). 

As mentioned in §0, the fragments formed from filamentary clouds tend to become prolate at the epoch 



of fragmentation. Therefore, as shown in §3.1 and §3.2, such a prolate fragment collapses to form an oblate 
core which thereafter evolves into a prolate core. This evolution is consistent with those of Figures § and ^. 



4. Time-Dependent Mass Infall Rate 



In this section, we discuss the evolution of the mass infall rate onto the center. 

Figure shows the evolution of the mass infall rate for model B. The mass infall rate is time-dependent. 
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in contrast to the prediction of similarity solutions (e.g., Shu 1977; Whitworth & Summers 1985). The mass 
infall rate takes its maximum at the very early stage and then monotonously decreases with time. This 
behavior is basically the same as that of Ogino, Tomisaka, Nakamura (1999)'s model with a small a, where 
a is the ratio of the gravitational force to the pressure force. Such evolution is explained as follows. During 
the runaway collapse phase, the central region tends to approach the Larson-Penston similarity solution 
(Larson 1969; Penston 1969) whose mass infall rate is temporally constant at Mi,p ^ Alc^/G. However, 
only the central small region converges onto the Larson-Penston solution by the epoch of core formation. 
In the outer infalling envelope, the pressure effect is important and the infall is appreciably retarded. Thus, 
after the central small region accretes onto the central point mass, the mass infall rate declines in time. 

Observed young stellar objects (YSOs) are often classified into several empirical evolutional stages from 
star-less dense cores to main-sequence stars (e.g., Andre, Ward-Thompson, & Barsony 1993). The youngest 
observed YSOs are called Class sources which are interpreted to be in a phase at which the envelope gas 
still has more mass than the central hydrostatic protostellar core. More evolved YSOs are called Class I 
sources. Recently, Bontemps et al. (1996) suggested that if the CO outflow rates are proportional to the 
mass infall rates onto the embedded young stars, the mass infall rates of Class sources are a factor of 
10 larger on average than those of Class I sources. Applying a simplified pressure-free model, Henriksen, 
Andre, & Bontemps (1997) proposed that the observed dispersion of the mass infall rates for Class sources 
is due to the time evolution of the mass infall rates. In our isothermal model, the evolution of the mass 
infall rate is qualitatively similar to that of Henriksen et al. (1997), although in our hydrodynamical model, 
the pressure force significantly retards the rapid change of the mass infall rate. As discussed by Ogino et 
al. (1999), the evolution of the mass infall rate also depends on the initial conditions of the clouds. We will 
discuss those effects in a separate paper. 



5. Implications to Observations 

Recent observations have revealed the detailed structures of the infalling envelopes around several 
embedded young stars (e.g., Hayashi et al. 1993; Ohashi et al. 1996; Saito et al. 1996; Momose et al. 1996). 
In this section, we discuss the physical structures of the disklike envelopes obtained by our calculations and 
make a comparison between our numerical results and observations. 

Figure ^ shows the radial profiles of p, v^., and in the equatorial plane during the accretion phase. 
Just after core formation, the density and velocity profiles are well approximated by p cx r^^, v^. ^ const.. 
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and ^ const, near the center. This indicates that the central high-density region converges well onto 
a similarity solution by the epoch of core formation (Matsumoto et al. 1997). As the collapse progresses, 
both the infalling disk and the rotationally-supported disk grow in size. In the infalling disk, the density, 
radial velocity, and azimuthal velocity have the radial distributions of p oc r^^ *^, Vr oc r^'''^, and oc r^^, 
respectively. The radial velocity is reproduced by a free-fall velocity of \/ 2GMq j r . The azimuthal velocity 
means that the specific angular momentum is almost constant in the infalling envelope. Such dependence 
seems to be qualitatively similar to that of a similarity solution of an infalling rotating disk (Saigo & 
Hanawa 1998; see also Terebey, Shu, Cassen 1984). In the rotationally-supported disk, the radial and 
azimuthal velocities are well approximated by ~ and v^p oc r"*^'^, respectively. At the boundary between 
the infalling disk and the rotationally-supported disk, a shock wave is formed and propagates outwardly. 

Momose et al. (1996) have examined the physical structure of the infalling disk around the L1551 IRS5. 
From a comparison between the observed postion-velocity diagram and their simple analytic model, they 
suggested that the velocity field in the disklike envelope is explained in terms of infall with slow rotation. 
Their simple model indicates that the infall and azimuthal velocities have radial dependences of oc r^*^'^ 
with 0.5 km s^^ at r — 700 AU and oc r^^ with 0.24 km s^^ at r = 700 AU, respectively. If we take 
po ~ 10^ cm~^ and T = 10 K, then the observed velocity field around the L1551 IRS5 is consistent with 



that of the model presented in §3^ at the stage at which i — icmo ~ 1.5(27rG'/9o) ""'^ — 3.8 x 10* yr, where 
t ~ ^corc denotes the time measured from the epoch of core formation. The mass infall rate is evaluated as 
M ^ 1 X 1O^^M0 yr^^. The central mass (a point mass plus its surroimding rotationally-supported disk) 
is then estimated as 0.6 Mq. 

The detailed comparisons between the numerical results and observations will reveal the formation and 
evolution processes of the infalling envelopes as well as the young stars. 
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Tabk 1. MODEL PARAMETERS 
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0.10 


15.5 


2.10 X 10^ 


1.52 X 10" 


-2 


1024 X 512 



Note. — P = the ratio of the centrifugal force to the 
pressure force; = the longitudinal wavelength; Mt = 
total cloud mass; Ar and Az = the grid spacings in the r- 
and 2:-directions; A^^ and A^^ = the number of grid points in 
the r- and ^;-directions. 




Fig. 1. — Density and velocity distributions in the r — z plane for model B at eight different stages. This 
model has the initial parameters of /3 = and Xz = 15.5ru. The gray scale indicates the density distribution 
in a logarithmic scale. The arrows show the velocity distributions which are normalized separately in each 
panel. Panels (a) through (h) show the cross sections at (a) t — tcoro = — 7.23iu (Pc = Po)i (b) —0.59 
(Pe = lOpo), (c) 0.00 (pe = lOVo), (d) 0.93 (Mo ~ 0.3Mt), (e) 2.25 (Mo ^ 0.5Mt), (f) 4.49 (Mq ~ 0.7Mt), 
(g) 8.28 (Mo ~ 0.86Mt), and (h) 10.88 (Mq ~ 0.9Mt), respectively, where t - Uore is the time measured 
from the epoch of core formation, pc is the central density, and Mq is the mass of the central point-like 
object. The first three panels represent the cross sections during the runaway collapse phase, while the 
others are during the accretion phase. Panel (i) denotes the enlargement of the central region of panel (h). 
If we adopt po = lO^cm"^ and Cg = 0.2km s~^, then the units of time, length, and mass are given as 
iu = (27rGpo)-°-5 = 2.5 x 10* yr, = Cst^ = 1.1 x 10^ AU, and = por^ = 7.6 x lO'^M©, respectively. 
The total cloud mass is equal to 1.9 x lO^rriu = 1.5M0. 



Fig. 2. — Same as Figure || but for model E. This model has the initial parameters of f3 — 0.05 and 
\z = 15.5ru. Panels (a) through (i) show the cross sections at (a) t — icorc = — 7.26iu {pc — po), (b) —0.59 
{pc - lOpo), (c) 0.00 {pe = lOVo), (d) 0.98 (Mq ~ 0.3Mt), (e) 4.71 {Mo ~ O.TMt), (f) 6.59 (Mq - 0.79Mt), 
(g) 8.25 (Mo ~ 0.82Mt), (h) 11.74 {Mq ~ 0.84A/t),and (i) 13.44 (Mq ~ 0.85A/t), respectively, where t-tcore 
is the time measured from the epoch of core formation, pc is the central density, and Mq is the mass of the 
central point-like object. If we adopt po = lO^cm"^ and Cs = 0.2km s~^, then the units of time, length, and 
mass are given as — (27rG'/9o)^"'^ = 2.5 x 10"* yr, = Cgt^ = 1.1 x 10^ AU, and m„ — por^ — 7.6 x lO^^M©, 
respectively. The total cloud mass is equal to 2.1 x lO^mu = I.6M0. 
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Fig. 3. — Density and velocity distributions of the central region for model E. Panels (a), (b), and (c) denote 
the cross sections at the same stages as panels (g), (h), and (i) of Figure ||, respectively, [(a) t — tcorc = 8.25iu 
(Mo =i 0.82A/t), (b) 11.74 (Mq ~ 0.84Aft), and (c) 13.44 (Mq ~ 0.85Mt), respectively, where t - tcore is the 
time measured from the epoch of core formation and A/q is the mass of the central point-like object.] 




Fig. 4. — Density and velocity distributions in the r — z plane for the prolate cloud at four different stages. 
This model has the initial parameters of (i?07-^o) = (1.66ru, 2.49ru). The gray scale indicates the density 
distribution in a logarithmic scale. Panels (a) through (d) show the cross sections at (a) t — tcorc = — 2.06tu 
(p, = po), (b) 0.00 (pc = lOVo), (c) 1.09 (A/o ^ 0.5A/t), (d) 2.64 [Mq ~ 0.7A//t), respectively, where t-Ure 
is the time measured from the epoch of core formation, pc is the central density, and A/q is the mass of the 
central point-like object. 
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Fig. 5. — Same as Figure ^ but for the oblate cloud. This model has the initial parameters of {Rq^Zq) = 
(2.40ru, 1.60ru). The gray scale indicates the density distribution in a logarithmic scale. Panels (a) through 
(d) show the cross sections at (a) t-t^ovc = -2.08tu [p^ = po), (b) 0.00 [pc = lO^'po), (c) 0.49 {Mq ~ 0.3Mt), 
(d) 1.12 (Mo ~ O.bMt), respectively, where t — tcorc is the time measured from the epoch of core formation, 
Pc is the central density, and Mq is the mass of the central point-like object. 
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Fig. 6. — Mass infall rate of model B as a function of time measured from the epoch of core formation 
{t — tcoie)- If we adopt po = lO^cm"^ and Cg = 0.2km s~^, then the units of time and mass infall rate are 
given as = (27rGpo)~°'^ = 2.5 x 10* yr and cl/G = 1.9 x 1O~^M0 yr~^, respectively. Open circles denote 
the stages at which Mq = O.lMt, 0.2Mt, • • •, 0.7Mt. For comparison, we show the cross sections of the cloud 
at the stages of Mq = 0.3Mt, 0.5 Mt, and 0.7 Mt. 
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Fig. 7. — Radial profiles of (a) density, (b) radial velocity, and (c) azimuthal velocity in the equatorial plane 
for model E. They are plotted at the stages at which (1) t — icoro = O.OOtu (the epoch of core formation), (2) 
1.56 (Afo — 0.4Aft), and (3) 13.43 (Afo — 0.85Aft), where t — tcorc is the time measured from the epoch of core 
formation and Mq is the mass of the central point-like object. If wc adopt po = lO^cm^'' and Cg = 0.2km 
s~^, then the units of time and length are given as tu = (27rG/9o) ""^ — 2.5 x 10^ yr and = Cst^ = 1.1 x 10^ 
AU, respectively. The total cloud mass is equal to 2.1 x lO^mu = I.QMq. 



